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ABSTRACT 
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\ 

^ I 

I In  a recent  paper  Charlock,  Herman  and  Zdunkowski  (1976)  disagree 

i I 

1 on  how  best  to  apply  the  Radiative  Transfer  Equation  (RTE)  in  the  con- 

I struction  of  Elsasser  type  radiation  tables.  The  two  proposed  approx- 

' imate  solutions  are  analyzed  and  compared  against  a quasi-exact  solu- 

I tion.  Infrared  fluxes  and  cooling  rates  are  calculated  for  part  of 

the  6.3  pm  water  vapor  band  for  two  model  atmospheres.  It  is  found 
I that  the  Zdunkowski  (Z)  approximation  yields  more  accurate  downward 

I fluxes,  while  the  Chari ock-Herman  (CH)  approximation,  in  general,  re- 

1 

suits  in  more  accurate  upward  fluxes.  For  the  two  model  atmospheres 
I studied  the  cooling  rates  for  the  Z approximation  are  usually  of 

better  quality  than  those  due  to  the  CH  solution,  unless  the  diver- 
I gence  of  the  net  flux  is  extremely  small. 
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I.  INTRODUCTION 


The  transfer  of  Infrared  Radiation  in  the  earth's  atmosphere  is 
of  great  interest  to  meteorologists.  Over  the  years  different  mathe- 
matical models  have  been  exploited  with  varying  degrees  of  success. 
All  of  these  are  special  cases  of  the  fundamental  Radiative  Transfer 
Equation  (RTE),  whose  solution  is  known  in  principle.  However,  there 
are  many  difficulties  in  their  practical  application  to  the  atmo- 
sphere. Hence,  various  approximations  are  proposed  to  simplify  the 
computation  of  atmospheric  radiative  fluxes  and  attendant  cooling 
rates. 

In  1966  Zdunkowski,  Barth  and  Lombardo  and  Zdunkowski  (1971), 
hereafter  referred  to  as  Z,  proposed  an  approximate  solution  to  the 
RTE  based  upon  the  original  work  of  Elsasser  and  Culbertson  (1960), 
hereafter  referred  to  as  EC.  In  1976,  Charlock  and  Herman,  hereafter 
referred  to  as  CH,  claim  to  have  obtained  the  "correct  expression" 
for  infrared  radiative  transfer  and  reject  Z's  approximation  as  in- 
correct. However,  Z (Charlock,  Herman,  Zdunkowski,  1976)  points  out 
that  their  expressions  are  merely  approximations  and  it  is  improper, 
a priori,  to  dismiss  one  approximation  in  favor  of  another. 

The  purpose  of  this  study  is  to  analyze  and  compare  the  atmo- 
spheric radiative  fluxes  and  cooling  rates  resulting  from  tli<t>two 
approximations, montioned  abov^  As  a norm  for  the  comparison,  a 
quasi-exact  solution  is  also  applied.  This  is  accomplished  by  using 
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the  EC  formulation  of  the  generalized  absorption  coefficient.  Calcu- 
lations are  carried  out  for  two  model  atmospheres  for  part  of  the 
6,3  micron  water  vapor  band. 

Initially,  it  was  planned  to  check  the  quality  of  the  two  approx- 
imations against  measurements,  but  it  is  believed  that  the  computation- 
al type  of  comparison  is  more  useful  and  definitive, 

""It  is  found  that  the  Z approximation  is  more  accurate  for  the 
downward  fluxes  and  the  CH  approximation  is  more  accurate  for  the  up- 
ward fluxes.  For  the  model  atmospheres  studied  the  cooling  rate  errors 
are  usually  smaller  for  the  Z approximation,  unless  the  divergence  of 
the  net  flux  is  extremely  small. 


II.  DATA  AND  ANALYSIS 


2.1  Symbols 

To  guarantee  a consistent  notation  it  is  useful  to  first  define 
the  following  symbols: 


B 


V 


7'  P T-T  P 
^ 'o’  0 


u = / Pydz’ 


T = / K^{u)du 


monochromatic  hemispheric  blacL  ivody  flux  at 
wavenumber  v;  a subscript  G refers  to  the  ground, 
height,  temperature,  pressure;  subscript  "o" 
refers  to  standard  values, 
absorber  density. 

absorber  mass  in  a vertical  column  of  air 
[gm  cm  ]. 

reduced  absorber  mass, 
absorption  coefficient. 

generalized  absorption  coefficient;  = I-^(Tq) 
optical  pathlength  (not  to  be  confused  with 
transmission  function), 
flux  transmission  function, 
downward  monochromatic  flux, 
upward  monochromatic  flux. 


t 


time. 
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I 


I 


f 


I 


2.2  The  Model  Equations 

For  a plane  parallel,  horizontally  homogeneous,  non-scattering 
atmosphere  an  exact  solution  to  the  RTE  exists.  Assuming  the 
absence  of  downward  radiation  at  the  top  of  the  atmosphere,  the  solu- 
tion to  the  RTE  for  downward  monochromatic  flux  (Figure  1)  at  the 

14 

reference  level,  u^,  is  given  by  (See  Appendix  I) 

^(ui)du^j|du  (1) 

The  monochromatic  upward  flux  at  the  reference  level  is  expressed  as 

' ^,G  I 

■ I K^(u,)du,^J  dG  (2) 

Note  that  the  absorption  coefficient,  K^,  is  a function  of  pressure 
and  temperature. 

2.3  The  Elsasser  Scheme 

In  order  to  facilitate  flux  calculations,  EC  replace  the  rapidly 
varying  absorption  coefficient,  K^,  by  an  empirically  determined, 
more  slowly  varying,  generalized  absorption  coefficient,  L^.  While 
varies  with  pressure  and  temperature,  is  assumed  to  depend  only 

on  temperature  (Figure  2).  The  temperature  dependence  is  given  by 


Figure  1.  Schematic  representation  of  flux  integration  limits 


7 


EC  (1960)  as 


log^gL 


- a 


(To-T) 


log 


10  T 


(3) 


The  quantities  L^,  and  a are  given  by  EC. 

In  view  of  the  EC  approximation,  the  pressure  broadening  effect 
on  the  absorption  coefficient  is  included  by  defining  the  reduced 
absorbing  mass,  i.e. , 


Using  Equations  (3,4)  the  argument  of  the  transmission  function  be- 
comes 


z “o 

/ [P(z’).T(z')]  p^(z')  dz'  = / [P(ui),T(u^)]  du^ 


Substituting  Equation  (5)  into  Equation  (1)  and  integrating  by  parts 
yields 


“o  . / “o  \ 


du 


8 


B^(T(u^))  - B^(T(u 


- / 
0 


*=0))Tp(/  L^(updu*^ 

(T(u*))  / % \ 


(6) 


for  the  downward  flux.  Similarly,  for  the  upward  flux  one  finds 


V “o  / % 


du 


"k 

^(updu^^du*  = B^(T(u;)) 
'‘^0  ' 


U*  d B^(T(u*)) 

^ t 

Uq  du 


(4  L^(updu;^du* 


(7) 


In  Equation  (7)  it  is  assumed  that  the  temperature  of  the  air  directly 
overlying  the  surface  of  the  earth  equals  the  surface  temperature. 
Transforming  Equations  (6,7)  in  temperature. 


★ 

fJ(u;)  = B^(T(u;))  - BjT(u*=0))  Tp(  /°  L^(updu;) 


d B^(T(u*)) 


((u**0) 


dT 


'f(/*  L^(uT)du;jdT  (8) 
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and 


F^(u  ) = B (T(u  ))  + / 

v'  o'  v'  ' o"  J(u*) 

'‘'o' 


u* 

(/. 


(9) 


Equations  (8,9)  contain  some  approximation  but  may  be  considered  quasi 
exact  flux  solutions.  For  simplicity,  they  will  be  referred  to  as  the 
exact  solution  hereafter. 

2.4  The  Chari ock-Herman  Approximation  (CH,  1976) 

Equations  (8,9)  refer  only  to  monochromatic  radiation.  Fluxes 
required  for  the  wavenumber  interval  (v-j , V2)  are  obtained  from 


and 


= / F*(u*)  dv 

J \)'  n' 


'’2  . 

/ F^(u*)  dv 

v'  0^ 


(.0) 


(11) 


The  radiation  table  approach  of  EC  requires  that  the  wavenumber  inte- 
gration be  carried  out  once  and  for  all  leaving  the  flux  for  a given 
atmospheric  path  a function  of  reduced  absorbing  mass  only. 

Inspection  of  Equation  (5)  shows  that  the  argument  of  the  flux 
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transmission  function  is  an  integral,  allowing  an  infinite  variety  of 
atmospheric  paths.  Clearly,  this  type  of  situation  cannot  be  accom- 
: modated  by  the  simple  radiation  table  approach  which  seeks  to  repre- 

sent fluxes  in  the  U*T-plane.  Therefore,  EC  introduce  a measured 
flux  transmission  function  defined  by  an  expression  of  the  form 

! 

i 

Tp  /*  = Tp  jL^(T(u*))[u*-u*] 


Euqation  (10)  assumes  that  is  evaluated  at  the  temperature  of  the 
elemental  emitting  layer  instead  of  evaluating  along  the  entire 
path.  This  assumption  is  appropriate  only  if  the  temperature  varia- 
tion along  the  atmospheric  path  is  very  small  (Sasamori,  1968). 

CH  (1976)  introduce  Equation  (12)  into  Equations  (8,9,10,11) 
and  obtain 

'’2  ( 

‘ I = B^(T(g*-0)  Tp(L^(T(u*=0))  up 

Vi  f 

d Bv(T(u*))  ) 

" ^ HT Tp(L  (T(u*))[u*-u*])dT  (13) 

T(u*=0)  ^ ° \ 

and 


2 \ * 

= / dv  Ib  (T(u*))  + / 
i I \)'  ' o''  *. 

V,  « T(Uq) 


T(U*)  d B^(T(u*)) 


Tp  (L^(T(u*(([u*-u*])dT 


2.5  The  Approximation  by  Zdunkowski , Barth  and  Lombardo  (1966) 


Z (Charlock,  Herman,  Zdunkowski,  1976)  looks  at  the  problem  from 

a different  point  of  view.  If  the  temperature  varies  considerably  a- 

long  the  absorption  path,  Equation  (12)  may  not  be  a sufficiently 

accurate  approach.  To  account  for  this  temperature  variation  a 

correction  term  Atc  should  be  added,  i.e., 
r,v 


'F  ( K ■ 

'll  ’ 


Tp(L^(T(u*))[u*-u*])  + 


^F,  (VU*.T(u*))  + ATp^^ 


Due  to  the  infinite  multiplied cy  of  absorption  paths,  it  seems  impos- 
sible to  incorporate  this  Atp  ^ in  the  radiation  tables.  Since  an 
approximation  like  Equation  (12)  must  be  made  in  the  radiation  table 
approach,  Z uses  a partial  derivative  in  the  fundamental  transfer 
equation,  e.g. , 
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1 


F^u*) 

V 0 


2 ( 3Tp 

= / dv  B^(T(u*))  ^ (L^(T(u*))[u;-u*]  du* 

vi  f 


(16) 


To  approximate  the  temperature  variation  expand  the  partial  derivative. 


. ^^F.v  ^ 

3u*  du*  du* 

Substituting  Equation  (17)  into  the  up  and  down  flux  equation  results 
in  an  approximate  solution,  used  by  Z,  ^ (1966),  of  the  form 


^2 

F^(u*)  = / dv j B^(T(u*))  - B^(T(u*=0))  Tp(L^(T(u*=0))u*) 

Vi  ( 


T(u*=0) 

d B^(T) 

dT 

- 

3 

B^(T)  - 

T{u*=0) 

V 

* . *1 


rF(Lv(T)[u;-u*: 

3T 


dT  (18) 


and 


T(U*)  d B^,(T) 


/ dv|B^(T(u;))  + / '.  ■ t,(L.(T)[u*.u:])  dT 


dT  ‘F'"v' 


T(U*) 

+ / ^ B (T)  ' ■ .V ^ dT 

Uu*)  ^ 


3tf(L^(T)[u*-u;]) 


3T 


(19) 


I 
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Comparing  CH  and  Z approximate  solutions  shows 


2 ( T(u* 

- / dvM  ° 

I r(u*= 


8T 


T(u*)  aXcCL  (T)[u*V])  ; 

B (T)  dT 

=0)  ^ 


(20) 


and 


= F, 


CH 


^2  / T(U*) 

«>-l,  ‘•jl,.;, 


3Tp(L^(T)[u*-u;]) 


3T 


dT 


(21) 


2.6  Numerical  Treatment 

The  reduced  absorbing  mass,  u*,  is  calculated  using  trapezoidal 
quadrature.  To  obtain  sufficient  accuracy,  u*  is  determined  at  every 
5 m for  the  nearest  100  meters  (m)  either  side  of  the  reference  level 
and  at  every  25  m for  the  remainder  of  the  atmosphere. 

For  each  of  the  solutions  (Equations  8,9,13,14,20,21)  to  the 
RTE,  the  monochromatic  upward  and  downward  flux  through  a reference 
level  is  calculated  again  using  the  trapezoidal  rule.  For  flux  calcu- 
lations the  nearest  100  m on  either  side  of  the  reference  level  is 
divided  into  25  m increments,  while  the  remainder  of  the  atmosphere 
is  divided  into  100  m intervals.  Thus,  for  example,  the  monochromatic 
down  flux  for  the  exact  case  (Equation  8)  is  calculated  in  the  follow- 
ing manner: 
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B^(T(up)  - B^(T(u*=0))  Tp 


L.)(u 


i+1 


N 

i=0 


d B 


■^Fd) 


v(i+1) 

dT 


^F(i+1)) 


(VrTi)j 


(22) 


In  Equations  (20,21)  the  partial  derivative  with  respect  to  tempera- 
ture is  determined,  as  appropriate,  in  one  of  two  ways.  At  the  sur- 
face of  the  earth,  the  top  of  the  atmosphere  and  at  the  top  of  an 
inversion  a one-sided  derivative  approximation  is  used.  For  example. 


■'i.i 


Ti 


(23) 


where  i and  i+1  are  10  m apart.  Note  that  the  reduced  absorber  mass 
is  held  constant  while  the  temperature  is  allowed  to  vary.  At  all 
other  levels  in  the  atmosphere  a centered  difference  scheme  of  the 
following  form  is  used: 


8 T ■ ■ ■ ' 


'i+1  - ^--1 


(24) 


The  slowly  varying  generalized  absorption  coefficient,  L , is 


-1 


tabulated  by  EC  for  40  cm  increments  of  the  wavenumber  spectrum  for 


j 
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a temperature  of  293°K.  For  other  temperatures,  is  found  from 
Equation  (3).  Since  and  d B^/dT  are  virtually  constant  over 
Av  = 40  cm“\  they  are  evaluated  at  the  center  of  the  increment. 
Fluxes  are  computed  as 


^(“o>  “ I 


N 

E F.(up  Av 


(25) 


The  radiative  cooling  rates  are  proportional  to  the  divergence 
of  the  net  flux,  f^,  i.e.. 


^ = d_  ^ A -J- 
dt  ~ PC^  dz  “ -r  Az 


(26) 


where  p and  Cp  represent  the  average  density  of  the  air  and  the 
specific  heat  at  constant  pressure,  respectively,  and 


F^  = F^z)  - F^z)  . 


(27) 


2.7  Properties  of  Approximations 

In  this  section  some  qualitative  statements  are  made  which 
describe  the  two  approximations.  As  previously  stated,  the  important 
difference  between  the  approximations  and  the  exact  solution 
(F^x(^p)  is  the  manner  In  which  the  argument  of  tp  ^ is  evaluated. 

First,  consider  the  downward  flux  in  the  Atmospheric  Model  I 
(Figure  3).  In  the  exact  case  one  obtains 


HEIGHT  (KM) 
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1 TEMPERATURE- DEW/FROST  POINT  (®K) 

I 

I 
I 


Figure  3. 


Model  atmospheres  I and  II. 
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,p(  ; L„(u;)dup  = 


Tr(L  (U*-U*)) 

F'  V 0 ' 


as  follows  from  the  mean  value  theorem.  On  the  other  hand,  the 
approximate  solution  evaluates  at  the  temperature  of  the  elemen- 
tal  emitting  layer,  i.e.,  Tp(L^(u*) [u*-u*]) . Since  the  temperature 
decreases  with  height,  it  is  evident  that  (Figure  2), 


Furthermore,  ip  ^ decreases  with  increasing  optical  pathlength  thus 


Tp(yu;-U*])  < Tp(L^(u^)[u;-U*]) 


*\r.  * ..*1 


By  comparing  corresponding  terms  in  the  exact  and  approximate  solu- 
tions for  downward  flux  one  finds  from  Equation  (30)  that 


B„(T(u*=0);  Tp  (/  \(T(up)(Jup  < B^(T(u*=0))  Tp(L^(T(u*=0) 


[u;-0]  ) 


T(uM  d B (T(u*))  “o  ^ ^ T(uJ) 

/ ° — TT </  LJT(dp)dur)  dT  < / ” 


T(u*=0) 


T(u*=0) 
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d B^(T(u*)) 
dT 


Tp(L^(T(u*))[u*-u*])  dT 


Comparing  Equations  (8,13)  it  is  evident  that 


(32) 


(33) 


for  an  atmosphere  in  which  temperature  decreases  with  height. 

Still  considering  downward  flux,  given  T2  > T-j  it  follows  (Figure 
2)  that 


Tp(L^(T2)[u;-u*])  - Tp(L^(TT)[u;-up)  <0 


(34) 


and,  therefore, 

3Tp(L  (T(u*))[u*-u*]) 

Hence,  from  Equations  (13,20)  it  is  seen  that 


(35) 


f5(u;)  > fJ,(u;)  (36) 

for  atmospheres  of  decreasing  temperature  with  height. 

Now  consider  briefly  the  upward  flux  for  Atmospheric  Model  I. 
Analogous  to  the  downward  flux,  it  follows  for  the  upward  flux,  from 
Figures  1,  2 and  3 that 
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4 


I 


1 


I 

I 

I 

I 

I 


i' 


L^{T(u*)) 


(37) 


and 


Tp(L^[u*-u;])  > Tp(L^(T(u*))[u*-u^]) 


(38) 


Therefore, 


T(U*) 

T(u;) 


d B^(T) 
dT 


u*  . T(U*) 

T,a  L^{T(up)du,) 


d B (T) 
v'  ' 

dT 


Tp(L^(T(u*))[u*-u;])  dT 


so  that  Equations  (9,14) 


(39) 


(40) 


Since  3Tp  ^/3T  < 0,  as  in  the  downflux  analysis,  from  Equations 
(14,21)  one  finds 


f 

CH 


(41) 


for  atmospheres  with  temperature  decreasing  with  height  throughout. 

Next,  consider  the  downflux  in  Atmospheric  Model  II  (Figure  3). 
Analyzing,  as  earlier,  the  exact  vs.  the  approximate  case  for  the 
region  above  the  top  of  the  inversion,  one  finds 

^EX^'^o)  ^ 
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On  the  other  hand,  for  the  upward  flux  at  reference  levels  below 
he  top  of  the  inversion,  analysis,  as  earlier  done,  shows 


L > L {T(u  )) 

V V 


Tp(L^[u*-u;])  < Tp(L^(T(u  ))[u*-u;]) 


Also,  since  now  ^ ^1  ’ follows  that 


Tp(gT2)[u*-u;])  >rF(^(Ti)[u*-u;^) 


so  that 


9Tp*L^(T(u*))[u*-up) 


Thus,  one  finds  from  Equations  (9,14) 


and  from  Equations  (14,21) 


1 4- 


III.  DISCUSSION  AND  RFSULTS 
3.1  Atmospheric  Models 

Two  model  atmospheres  are  considered  in  this  study.  The  first 
(i.e.,  I-Figure  3)  is  the  NACA  Standard  Atmosphere  with  an  assumed 
dewpoint  distribution  typical  of  a cloud  free  atmosphere.  The  upper 
level  moisture  profile  from  7.2  to  12  KM  uses  a representative  frost 
point  temperature  sounding  (Craig,  1965).  The  upper  and  lower  mois- 
ture profiles  are  matched  in  such  a way  as  to  result  in  a continuous 

function  of  water  vapor  density.  The  total  reduced  absorber  mass  of 

_2 

the  atmosphere  is  approximately  1.2  gm  cm  . The  second  model  atmo- 
sphere (i.e.,  II-Figure  3)  is  similar  to  the  first,  but  a moist  inver 
sion  is  introduced  in  the  lowest  500  meters. 
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3.2  Flux  Transmission  Function 


■1 


The  spectral  region  investigated  stretches  from  1100-1420  cm 
and  contains  part  of  the  vibrational-rotational  water  vapor  band  and 
the  absorption  window.  In  general,  the  flux  transmission  function 
which  consists  of  two  parts  is  given  by 


^F, total  * ^F,cont.  ^F.band 


(50) 


Tp  is  evaluated  in  terms  of  the  generalized  absorption  coefficient, 
L^,  while  the  flux  transmission  function  for  the  continuum  is  expressed 
in  terms  of  the  continuous  absorption  coefficient,  k , i.e.. 
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Tp  IJ*)  ” 2E.  (k  U*)  = 

F.cont  V V 


-1 .66k\,  u 


(51) 


E.  is  the  exponential  integral  of  order  3.  It  should  be  noted  that 
3 

is  multiplied  by  a factor  of  ten,  as  implied  by  EC  (1960). 

The  flux  transmission  function  for  6. 3 u band  is  approximated  by 
a piecewise  cubic  spline  interpolant,  i.e.. 


TF^y(x)  = C^.  + Ax(C2^  + Ax(C3.  + AxC^.))  (52) 


X is  the  optical  pathlength  whereas  €2^,  C^.  and  are  inter- 
polant coefficients  determined  for  each  interval,  i,  of  data  points. 
The  quality  of  the  fit  to  EC  data  is  quite  good,  as  verified  by 
Table  1. 

3.3  Infrared  Fluxes 


Upward,  downward  and  net  fluxes  for  Model  Atmosphere  I (Figure  3) 
are  given  for  selected  reference  levels  for  the  spectral  region 
1100-1420  cm’^  in  Tables  2,  3 and  4.  Also  shown  in  the  Tables  is  the 
error,  E,  of  the  approximate  as  measured  against  the  exact  solution, 
i.e. , 


E(%) 


^approx^%^  ~ ^EX 


(53) 


Table  2 indicates  the  upward  flux  is  generally  more  accurate  for 
the  CH  than  for  the  Z approximate  solution.  This  is  in  agreement  with 
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Equations  (40,41).  Table  3 verifies,  on  the  other  hand,  that  the  Z 
approximation  yields  downward  fluxes  generally  more  accurate  than  the 
CH  solution.  This,  again,  satisfies  the  inequalities  of  Equations 
(33,36). 

Table  4 shows  that  the  net  flux  error  for  the  Z approximation 
is  less  than  6%  for  the  selected  levels  chosen  while  the  CH  approxima- 
tion yields  maximum  errors  of  almost  12%.  The  CH  approximation  is, 
however,  generally  more  accurate  in  the  upper  levels  where  the  down- 
ward flux  is  small. 

For  the  Model  II  Atmosphere,  as  seen  in  Table  5,  the  CH  approxi- 
mation, again,  leads  to  smaller  upward  flux  errors,  yet  the  maximum 
Z error  is  less  than  3.5%.  The  relatively  small  Z errors  result  from 
more  accurate  fluxes  in  the  inversion  region  as  qualitatively  expected 
from  Equations  (48,49). 

Table  6 shows  that,  as  in  the  Model  I Atmosphere,  the  Z solution 
yields  much  more  accurate  downf luxes  for  Model  II.  The  CH  approxi- 
mation produces  errors  of  greater  than  10%  at  all  levels  shown. 

The  net  flux  calculation  for  Model  II  are  given  in  Table  7. 

Again,  it  is  evident  that  the  CH  approximation's  accuracy  is  best  in 
the  upper  levels,  while  the  Z solution  is  best  in  the  low  levels. 

Yet,  over  all  levels  the  Z method  yields  a maximum  error  of  less  than 
3.5%,  whereas  the  CH  method  yields  a maximum  error  of  less  than  15.5%. 
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3 


31 
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3,4  Cooling  Rates 

Cooling  rates  are  computed  in  degrees  per  day  for  selected  layers 
of  100  meter  thickness.  Table  8 indicates  the  cooling  rates  for 
layers  centered  at  the  heights  shown  for  Model  I.  The  Z cooling  rates 
are  more  accurate,  except  at  the  top  of  the  atmosphere  where  the  heat- 
ing rate  is  vanishingly  small.  Table  9 shows  that  the  accuracy  of 
the  Z cooling  rates  for  Model  II  is  better  than  those  of  CH,  except 
at  the  top  and  bottom  of  the  atmosphere,  where  the  cooling  rates  are 
extremely  small.  This  suggests  that  the  Z approximation  yields  more 
accurate  flux  divergence  values  except  in  regions  where  the  divergence 
is  vanishingly  small,  such  as  at  the  top  of  the  atmosphere. 


i 

I 
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IV.  CONCLUSION 


It  is  apparent  that  the  CH  approximation  will  yield  more  accurate 
upward  fluxes  in  atmospheres  of  decreasing  temperature  with  height. 
Equally  apparent  is  the  fact  that  the  Z approximation  produces  more 
accurate  downf luxes  in  such  atmospheres. 

The  net  flux  through  reference  levels  computed  by  the  Z approxi- 
mate solution  seems  to  give  a smaller  error  than  that  by  CH.  However, 
for  the  upper  levels  of  an  atmosphere  the  CH  approximation,  in  general, 
yields  more  accurate  net  fluxes.  If  flux  divergence  or  cooling  rates 
are  desired,  it  appears  that  the  Z approximation  gives  better  accur- 
acy except  in  regions  where  the  cooling  rates  are  vanishingly  small 
and  likely  have  no  real  significance. 


I 


L 

APPENDIX  A 

RADIATIVE  TRANSFER  EQUATION 

To  guarantee  a consistent  notation,  it  is  useful  to  first  define 

the  following  symbols: 

4. 

I^  downward  intensity  at  wavenumber  v. 

I^  upward  intensity  at  wavenumber  v. 

u COS0,  9 is  the  zenith  angle  of  the  radiation 

4)  azimuthal  angle  of  the  radiation 
4- 

downward  monochromatic  source  function 
upward  monochromatic  source  function. 

A convenient  and  sufficient  starting  point  is  Schwartzschild's 
Equation  for  downward  radiation, 

dI^(T,u,4i)  , , 

u dr ^ J^(t,u,(}))  (1) 

and  for  upward  radiation, 

dI^,{T,u,4')  . + 

u = + Iy(T,y,(|))  - J^(t,u,i1»)  (2) 

A linear  differential  equation  of  the  form, 

^ + P(x)  y(x)  + Q(x)  = 0 (3) 

; has  the  following  solution 

I 

i 
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P(2)clz 


/ P(z)dz 


y(x)  = e Ly(x-)  - / Q(t)  e ° dt].  (4) 


Thus,  from  Equations  (1,2)  one  finds 


I^(T,y,(i>)  = I^(t^,u,(1))  exp(-(T^-T)/u)  + 


' I 

/ jUt,Us<P)  exp(-lt-T)/u)  ^ 


(5) 


for  upward  intensity  and 


I^(T,y,(t>)  = 1^(0, y,()))  exp(-T/y)  + 


dt 


/ J^(t,y,4>)  exp(-u-t)/y)  — 


(6) 


for  downward  intensity  where  from  Figure  4, 


T = / K^{ui)du^  i optical  pathlength  for  the  reference 

level  (7) 


U 

= I K^(u^)du^  H optical  path  for  an  arbitrary  level  (8) 


^Iv  " / ~ total  optical  path  for  the  atmosphere 


(9) 


V 
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See  Figures  1,4.  Under  the  Black  Body  assumption  the  infrared  inten- 
sity is  defined  by  the  Plank  function 

I^(p,(j,j  = B^/ti  (10) 

Therefore,  the  radiation  fluxes,  as  obtained  from  the  intensity 
by  angular  integration, 

2ti  1 1 

- ! i ydy  = 2 / B^udy  (11) 

0 0 0 

Transforming  Equations  (5,6)  using  Equations  (10,11)  yields  the 
monochromatic  flux  equations.  Assuming  the  radiation  at  the  top  of 
the  atmosphere  is  zero,  the  monochromatic  downward  flux  at  a reference 
level  is 

T 1 

F^t  ,)  = 2 /^  / B (t  ) exp  {-  1 (x  - t )}  ydy  ^ (12) 

v v * V V y V V y 

0 0 

The  monochromatic  upward  flux  at  the  reference  level  is 

2 ^,G  I i (■'iv  - 

■^Iv  1 

+ 2 / B^(t^)  / exp  {-  1 (t^  - T^)}  ydy  ^ (13) 

TO 

V 

Differentiating  Equation  (8)  with  respect  to  the  upper  limit  yields 

dt  = K (u)  du 

V V 


(14) 


Figure  4.  Schematic  representation  of  optical  path  for  flux  integration 
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Also, 


Uq  y Uq 

^ " / K^(u^)du^  - / K^(u^)du^  = / K (u^)du^  (15) 

0 0 u 


Substituting  Equations  (14,15)  into  Equations  (12,13)  gives,  for 
the  downward  flux  at  a reference  level. 


o 1 ^0 

= 2 / B^(u)  / exp  { - 1 / K^(u^)du^)  dpK  (u)du 

0 0 ^ u 


(16) 


and  for  the  upward  flux. 


f 1 , U 

" ^^v,G  ^ -jj  / K^(Ui)du^}  duK^(u)du 


U 1 , u 

2 / B (u)  / exp  {-  ^ ! K (u,  )du, } dpK  (u)du 

0 ^ 

0 0 


(17) 


It  is  advantageous  at  this  point  to  eliminate  the  integration 
overy.  Recall 

1/y  = sec  0 i C then  d^  = sec  9 tanO  d0  and  dy  = - sin  0 d0. 


Let 


X = / K^(u^)du^ 


(18) 


Then 


40 


1 , '^0  ti/2 

2 / exp  { - - / K^(u^)du^l  dy  = 2 / exp  {- 


/ K^(u^)du^}  sine  de  = 2 / exp(-xC)  ^ • (19) 


Introducing  the  exponential  integral, 


E.  (x)  = / exp(-x  ) ^ 


produces 


1 , “0 

2 / exp  ^ / K^(u^ )du^}  ydu  = 


2E.  ( / K (uJdui)  (21) 
^2  u ^ ‘ ‘ 


The  derivative  of  the  exponential  integral  is  defined  as 


dE.  (x) 
n 


= - E,  (x) 


Thus  it  is  evident  that 


- E.  (x) 
2 


dE.  (x) 
du  dx 


= I ^(iJi)dui)  = - 


Substituting  Equation  (24)  into  Equation  (23)  and  rearranging  yields 


Now  as  is  customary  define  the  flux  transmission  function 
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dE,  (x) 


(25) 


V,v  = 2E.^(x) 


(26) 


From  Equations  (21,25,26)  and  Equations  (16,17)  one  finds  for 
the  downward  flux  at  a reference  level 

u u 

^ fV.v<  / ’‘'“l » '1''  <”) 


Similarly,  for  the  upward  direction  at  a reference  level. 
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APPENDIX  B 

LISTINGS  OF  COMPUTER  PROGRAMS 
Spline  Interpolant  Program 


«• 

;>• 

o» 

7» 

t>* 

9» 

10* 

ll» 

12* 

13« 

IH. 

lb* 

16* 

17* 

16* 

IV* 

23* 

21* 

22* 


Ui.ltNSiON  XI  IbOl  iC  (9,50)  .Xl(  197)  >7(197)  tNPI  1)  imI  (SI  ,A2(5I  >AJ(b) 
(jrtT^  (Ai  1 1 ) , i = i,bi/>  cubic  splike  Interpolant  •/ 

EimTa  (a2(  I I , 1 = 1 .bl/'  X •/ 

uATA  (AJ(I).I=l,b)/'  FX  </ 

NPd  1 = 197 
REAo  (b.SOOlM 
bOO  FUR|-.aT(I2) 

RE AC ( 5. be  1 1 (XI ( 1 1 ,C(1, 1 1 > 1 = 1 ,N1 1 

C(2,1)=C(1.1)*(2.*XI(1|-X 1(2) -XI(3|)/((XI(I)-XI(2)|»(XI(1)-X Kill) 
•♦C(l,2|.(Xi(l)-x;(bl)/((XI(2)-xi(l||»(Xl|2)-XI(3)l) 
•*C(l,3).(Xl(ll-Xl(2))/((XI(3)-XI(ll).(XI(3|-XI(2||| 

C(2.r,l  |=c  ( 1 ,N1-2I  * (Xt  (hll-XI  (1(1-1  I I/I  (XI  (N1-2I-XI  (Ml-I)  I*  (X|  (Nl-2) 
l-XI(M)  ) i+ci  I >NI-1I»(XI1N1  l-XI  (01-2)  )/(  (XI  (N1-1)~XI  (Nl-2)  1* 

1(XI  (N1-1)_XI  (Nil  ) I ♦Cd.Nl)  » (2..XI  (la  l-xi  (N1-2I-XI(h1-1)  )/( 

1 (XI  (Nli-xi  (Nl-£)  I»(X1  (M)-XI  (Nl-ll  I ) 
n=U1-1 

CALL  SPLINEIN.Xl.C) 

CALL  CaLCCF(N,x1,CI 
bOl  F0HR«Ti2F5.2) 

X = AI ( 1 1 
uU  10  J=1.197 
X3( j|=X 


23* 

29* 

2b* 

26* 

27* 

2t* 

29* 

30* 

31* 

32* 

33* 

39* 

35* 


600 

10 

601 


FX=PcUolC(X,N,xl,C) 

7 ( J) =FX 

«R*TE(b,bOOI  J.X.FX 
FOHf at ( I 5,F10. 5,220.91 
X=a*,025 
UO  bOl  lrl,l.l 

‘•rite  (7, 601)  C(2,II  ,c(3,I),C(9,i) 
FOk|.,aT  (3Llb.9) 

CAlL  lL>f'L0T(12.,I2.) 

CALL  £2PlOT(X3,T,1,UP,1,1,a1,A2,A3) 

call  FINI 

STOP 

END 


Spline  Program  (continued) 


!• 

2« 

w* 

4* 

b* 

O* 

?♦ 

d* 

10* 

!!• 

12* 

U* 

iH« 

ib« 

ib« 

17* 

ib« 

iV* 


SUuPOUTInE  S>*LlNe(N.Xl*C) 

•jI«LiJS1Cn  XI  (bO)  .C(4,50)  »Di5oj  ,DIAG(50) 
LAff  Oi AO(l  I *01 L)/l • *0,  / 

NPXtNti 
UO  10 

10  ui (w » = (c ( I *M)-c  1 1 •^.-l n /o(M) 

uO  20 

C(2,Wl=J.»(Gltl|.&U&(M+l)*D(Mtl)»DitG(K.)  1 
2C  U1mv>('^)=2.«(LI(I-<I»U(H*1)  ) 

Uu  iO  1-1=2  >N 
G=-i,lNtl  )/D1*G(M-1) 

JO  C ( 2 « ,'0  ) =C  ( 2 • M 1 40  «C  ( 2 t I ) 

NJ=I,P1 

DO  40  .-.=2>U 

NJ=,.o-l 

40  C (2t,4j)  = (C(2»Nj)-D(hj)  •C(2.MJ'H  ) I/OIAGINJ) 
KCTGi^N 
El«0 


!• 

2> 

J« 

44 

5» 

6* 

7» 

a*  10 

9» 

10* 


1« 

2« 

i* 

•*9 

b9 

t9 

10 

19 

09 

i0« 

19 

11* 

12« 

20 

14* 

ib* 

JO 

it9 

179 

bUURourii.E  c«LccF(r.,»i,c) 

OlMEtiSION  XI<S0liC(4,S0l 
00  10  1=1|N 

ox=xi(i4ii-ii(ii 

i)IVCFl  = (C(l.l*l)-C(l.I))/Dx 

oivorJ=C(2. 1 )4C(2i  »♦!  )-2.»oiv:,fi 

C ( J . I ) = ( 0 1 VCM  -C  1 2 < I ) -u  I VDF  -5 ) /DX 
C(4,I)=OlVOFJ/UX/OX 

KETukN 

ENO 


FU!,crlull  PCUoICiXOAR.N.XI.C) 

OlyENSlUN  XI (50) >0(4,50  I 
OAIA  1/1/ 

OX=XuAH-XI ( I ) 

IF  (OX)  10.0020 
IF  (I. £0,1)  GO  TO  Jo 
1=1-1 

l)X=XUAi)-XI  ( 1 ) 

IF  (JXI  10. 00, JO 

l = l4l 

UX=UJX 

IF  (2.E0.N)  GO  TO  Jo 
00X=X8aH-X((U1) 

IF  (oOx)  00,19.19 

PCUbIC=C(l.n4ljX«(C(2.1l40X»(C(3.Il40X*C(4,l))  I 

HE  TORN 

Ei)0 


BKI  AVAIMSIE  COPY 
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Flux  Program  for  Model  II 


1* 

OlMLt.SlCN  UtROO  1 ITE  (900  ) .XT  (50  ) ,C  (4,501  ,GLZ(900  ) 

2* 

U:50 

3* 

«C AO (5. 500)  1C(2,I),C(3,II,C(4,iI,c(1,I),XI(I),T=1>N) 

4* 

500 

FOR  <AT (5C15.9) 

b* 

50? 

REA0(5,‘'Oi  . E (10=990 )ZRL.V,GL0> CON 

to* 

501 

TOHUaT (FG.n,F5.0,Fb.2>El2.fl) 

7* 

02=12000. 0-ZRL 

8* 

>'C=o2/100. 

9* 

xo=.-'3-l 

10* 

“u=2nL/ino. 

11* 

MU=r'U- 1 

l^* 

202=12000. 

Xb* 

22=202 

IH* 

L = 1 

lb* 

U(1)=0. 

lb* 

T2( 1 )=T (202) 

i7* 

CL2(1):GL(T2(1).V.&L0) 

18* 

c compute:,  u*  at  each  level 

19* 

IF(HC.LT.O)  uo=o. 

^0* 

IF(HO.LT.O)  Td=T(2C2) 

21* 

IF(HO.LT.O)  go  TO  8 

22* 

IF(MO.E«.0)  go  to  9 

23* 

00  1 I = l>Mr) 

244 

DO  2 J=l,4 

2b* 

L=L*1 

2b* 

21=22-25. 

27* 

T2(L )=T(21 ) 

26* 

GL2(L)=GL(T2(L)  ,V,gLO) 

29* 

UU:.5*  (01(0(22)  • (P(22)  21013.  )•  ( (293./T(22)  ) ••  ( .5)  ) ♦RH0(21  ) ' 

30* 

11013. )• ( (203.2T (21 ) ) ••( .5) ) ) • ( 22-2 1 ) • 1 00 . 

3t* 

U1L)=U(L-1)*UU 

32* 

22  = 21 

33* 

2 

COUTInUE 

34* 

1 

tONtPlUF 

3b* 

c coMpuits  u*  fOP  iiiIehval  ion  m above  reference  level  to  loo  m 

3o* 

22  = 2RL‘  mo. 

37* 

9 

00  3 1=1*5 

I 
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id* 

If 

O 

O 

is* 

L;L.  1 

HO* 

HI* 

TZ(Ll:T(71) 

H2* 

GLZ  IL )=CL ( IZII ) . V.oLOl 

Hi* 

UU^.b  - (FHOIZZ  ♦RtiO(Zl  ) • (P(2l ) / 

HH* 

1 1013.  )•  1 121  1 ) ••  1 .b)  ) 1 • CZ2-21  ) • 100. 

Hb* 

U(L)=U(L-H  ^UU 

H6* 

2?=21 

H7* 

“ COMTIMUC 

He* 

3 COliTIMUF 

49* 

UOFUIL  » 

50* 

TOFTZa  1 

51* 

IF(liU.LT.O)  GO  TO  11 

52* 

8 DO  7 1=1 .b 

5i* 

DOlO  JFl*** 

5*** 

LiL*l 

55* 

21=22-5. 

56* 

T2(Ct:T(21 1 

57* 

GL2 (L )=GL IT2(L) . V.uUO) 

56* 

UU=.5» (PH0(22> • (PI22)/1013.)*( (203./T (22) ) •• ( .51 ) ♦RHO(2l ) • (P (Zl ) / 

59* 

IIOIO.I^IIZTS./KZIII'M.S)))  *(22-211  *100. 

60* 

l;(L)=u(l-i  )*uu 

bl* 

22  = 21 

b^* 

10  CCMTIMIF 

bi* 

7 CONTIIlUE 

04* 

c 

computes  U*  poo  remaining  UEVELS  down  to  the  surface 

bS* 

22:ZRL-inO. 

6b* 

IF (MU.EO.O)  go  to  11 

b7* 

00  b I=1>MI) 

bd* 

00  0 J=l  -I* 

69* 

LrCH 

70* 

21=22-25. 

71* 

T2(L)=T1Z1) 

72* 

GL2(L):cL(T2(L) .V.OLO) 

7i* 

UU  = .5* (BHO( 22 )*(P(Zc)/1013.)*( (2P3./T (22) ) •* ( .5) ) ♦RHOIZl ) • (P(Z1 ) 2 

7H* 

11013.)*((2a3./T(Zl))**(.5)) )• (22-21 ) *100. 

75* 

0(L)iU(L-l)»UU 

7b* 

Z2rZl 

77* 

b CO(tTINUE 

76* 

5 CONTINUE 

79* 

c 

COMMUTES  U*  AT  THE  SURFACE(UG) 

60* 

UG=U(L) 

81* 

c 

COMMUTES  the  down  Flux  at  ThE  reference  level  for  ch  ano  z methods 

62* 

11  1FLAG=1 

6i* 

J=1 

6H« 

IFiMO.LT.O)  0FTC=0. 

85* 

IFIHD.LT.OI  DFTZ=0. 

66* 

if(mo.lt.O)  go  To  10** 

67* 

TS=T2( 1 ) 

86* 

ULr&L (TS> V.OLO) .ALOOIOIUO) 

89* 

If (V.GT.I2A0.)  COn:0, 

90* 

THC0N=EXP(-10.(.*CnN*U0) 

91* 

lHX=rcUP  IC  (UL  .11)  •TXCON 

92* 

Of  U=(!{  TA  (TO.VI-OEIAI  T2(  1 ) ,V)  ‘Tnx 

9i* 

GLUGLITS. V.OLO) 

9h* 

i 0=UtTA( To.V) 

9b* 

eir:ir TA i 1 1 1 . vi 

90* 

white  16,^0't  lUL  . THX  >PFt! 

97» 

WRITE  IU0<0I-1  tHO<Ul  >DFB 

9a* 

620 

F0i<"AT(5)(£  1^  .fli  1 

99* 

f04 

FOH  ■At ( ir  . XIElfa.S) ) 

100* 

OF  lCr-0. 

101* 

bfZ-0. 

102* 

r?xt?  1 1 » 

103» 

IFIfO.Nf.O)  GO  TO  112 

10*4* 

Tm  IrTHX 

10b* 

0B2:DH(T2.V) 

1 Oo  • 

B2li)tTA|T2,v) 

107* 

DTH2:oTP( I2.uO>N<2O2<V.0LO.2FL.c0n) 

loa* 

GO  TO  103 

109* 

112 

00  100  ‘'rl.'-o 

110* 

JrJ*4 

111* 

201=202-100. 

112* 

Tl=T2lJI 

113* 

U1=U0-'J(  J) 

114* 

IFIP.GT.I I GO  To  105 

115* 

TRIlrTR* 

11b* 

0TR2=CTPI T?,uO.N'202. V.GL0.2nL.CON) 

117* 

002=DB(T2.V) 

lie* 

32:bCTA(T2,vi 

119* 

105 

OLl=GL2(J)*AL0G10tUl) 

120* 

TRCCIf;I=f  XPl-lG.G^COf/'Ul  > 

121* 

TRI  = TRCC'II»PCUalC(ljLl.N) 

122* 

Bl=nETA<Tl,v) 

12J« 

OBl=Of)(Tl  .iM 

124* 

OF  = .5*  irn2»TRi i*ori»tri i • i ti-T2) 

I2b* 

DFIC:0FIC*0F 

126* 

107 

DTRl:orR(Ti,ui.n.20l.V,GL0.2RL,CONl 

127* 

DDF2P.5* 1D’*DTR2*r1«0TH1 ) ♦ 1 T1-T2) 

124* 

DF2=0FZ  ‘Onpz 

129* 

202=201 

UO* 

WHITE  ((..FOR)  J.OTRl  . TRI  .OOFZ.OF.zm  .01 .6L2(  Jl  .UL1.COM 

lil* 

690 

FORPAT  (IX.n.eiElR.Ol.lX.eiZ.B) 

132* 

triutpi 

133* 

OB2=D01 

134* 

B2=B1 

13b* 

OTR2=OTH1 

136* 

100 

T2  = T1 

137* 

103 

00  101  P=J.5 

13H* 

JZJtV 

139* 

201=Z02-?0. 

140* 

Tl=TZtJ) 

141* 

01=00-01 Jl 

142* 

106 

IFlOl.Lf.O.I  GO  To  102 

143* 

ol:=glzij) >ALOGiaioi i 

144* 

00  TO  IPO 

I4b* 

102 

0E1=-3.F 

146* 

109 

lRCOHI=f  XP(-16.6*C0M«01 1 

147* 

TR  UIHCP'iI  .CCOUIC  lULl  .Ml 

t4C* 

Otll=C'HtT  1 .VI 

149* 

Op=.b. (rn2*Tni I .OMj.TRl 1 • ( T 1-T2) 

IbO* 

Ibi* 

Of  ICrOFIC.'-F 

47 


IbJ* 
lb4  • 
IbS* 
Ibb* 
ib'* 
iba* 
lb4* 
loll* 
Ibl  • 
16^* 
lb3* 

164* 

lob* 

lob* 

lb7* 

lbb» 

lb9* 

170* 

171* 

172* 

17^« 

174* 

I7b* 

17b* 

177* 

17b* 

17'^t 

IbO* 

lei* 

ld2* 

Idi* 

164* 

las* 

16b* 
167* 
16b* 
189* 
19C* 
191* 
19^* 
19J* 
l<tHm 
I9b* 
19b* 
197* 
I9d* 
19'>* 
200* 
201* 
202* 
20^* 
2u4  • 

2ob* 

2i/b* 

207* 

2(it* 


f3l=^)r.lA  / Tj  ,v> 

lOrt  UlWirDTW  (Tl  .Ul  »rj*20l  • V.rA0*29t  • l^LAO' 

0DF2r  .5*  <07»t  »D7*» I ) • I Tl-T?) 

*0nf^ 

.«WnE(b.b9»IJ*0TRl.TOr#0**'P2»nr,/01*Ul.GL;^(J>»ULWC0N 

IF  < J.LT.Ub'>.Or<.j.r,l  .477)  GO  TO  UO 

TL  (b.bo;^)  J.202.  I2.COf/»  fwj  I »OT»2»f>f' 

I Tf.  (6#bO?)  J»OF  IC  •P^TP'OhI  #02*01  *0^2 
I Tt (b.bO?)^*20l » 1 1 •ULl 'OL/^ JJ  »TRt .OTRl 
I Tf.  fbf^O  ) #GL2{  J-4  ) 

b02  FORMAT  ( 1 X . I )«b(  1 < *L  16#  AH 
603  FORMAT  ( 1 X . 13*  3(Elb. AM 
no  TrIIsTRI 
e2=Bl 
0B2r0Bl 

0TR2=0TrI 

T2  = TI 

101  202=201 

OFTc=r>ro-oFic 
OF  T2rOFR-OFlC-OF2 

c COMPUTES  The  up  flux  aT  the  FEFEPENCE  level  for  The  cm  ano  z cases 

104  UF(C=0. 

UFZ=0. 

IFLAGS2 

UFb=OETA(TO  *VI 
IF (MU.LT.O)  go  to  20S 

c computes  thf  uppLux  for  The  layer  from  The  Reference  level  to  loo  m SEUdw 
00  200  Mri.S 

J=J*4 

201=202-20. 

Tl=T2< J) 

Ul=U< J)-00 
IF(M.GT.I>  go  to  206 
TRim . 

T2=T2(J-4) 

OB2=OB^T2.V> 

02=t5ETA(T2.V) 

OTR2=0. 

206  ULl=r.L2(J)*ALOGlO(Ul) 

TflC0Ml=EXPl-l6.6*C0N*Ul I 
thi=thconi*pcu0ic(uli.n) 

OenOBlTl  »V) 

UFr  .5*  <pP2*T«n  *091  * TR]  > • ( T 1-T2 ) 

UF 1C=UF1C*UF 
Dl=UETA(ri, V) 

207  OTRl=OTnn  1 »Ul  •N*20l  *V.r,LO*ZRL*cON» 

UUF2-.S* (R?*DTR2*ra  *0 TH I i • ( Tl -T2 > 

UF2=UFZ*UUF2 

Jj2=J-9 

wRi T£  <6»6Pa) J»0TRI •TRI .UUF2.UF.20l »Ul *0t2( J1  *ULl»C0N 
IF  { J,LT.47r|,0R.j.r,T.4«5)  GO  TO  2ll 
*R  I Tf  (6.60^1  J.  202.  TRI  I »0TR2.UF 

1 IK (6.607) J.UF IC>0U2*0P I .n2»nl *UF2 
NfM  TE  (6,607)0.201  . T I #((L  I •OC/T  J»  . TRf  .OTRl 
*iRl  11  (6,60  3)  JJ2*U(J)  »U(  J-4  ) »GL2(U-<«) 

211  ThM=TH| 

0u2=E(il 


w miUBtE  COPY 


209* 

H2I01 

210* 

0TK2r0TRl 

2ll* 

T2FT1 

212* 

200 

iOi-ZOl 

2i  i* 

C 

coMPuits  The  up  elux  fob  The.  lxyer  from  loo  m i 

2 U* 

c 

surface 

2lb* 

IF (MU.EO.OI  GO  TO  205 

2 16* 

Co  201  MFl.MU 

217* 

201=202-100. 

23ft* 

JZJtU 

219* 

Tl=T2tJ) 

220* 

U1FU( J»-UO 

221* 

UL1=GL2(J)‘AL0610IU1) 

222* 

0(il=DO<TI  .VI 

225* 

THCUfll=FXP(-i6.6*C0M*Ul  1 

22«4* 

THI=TRCTll«PCUOtC(ULl  .Ml 

225* 

209 

UF=.S* ICP2»TBI I»0R1»TRI 1 ‘(Tl-T?) 

220* 

UF  1C=UF1C*IJF 

227* 

Bl=bETA(Tl.V) 

22rt* 

210 

OIHirOTR ITl .U1 .N./Dl •V.0L0.2HL.C0NI 

229* 

UUF2=.5« (B2»0TK2*R1*0TR1 ) •(T1-T2I 

230* 

UF2=LF2.UUFi 

2il* 

■ H ITE  J'CTRI  . TR  t .UUF2.UF.  201.01  .6L2 

232* 

Th11=TR1 

2ii* 

0B2=DB1 

25*** 

92=U1 

23b* 

OTH2:oTR1 

23o* 

T2  = T1 

237* 

201 

202=201 

236- 

205 

uFrcruFo‘UFrc 

239* 

UFT2=oF.X*UF  ic*  UF2 

2hU* 

c 

computes  the  oowii  flux  at  tmc  reference  level 

2>4l* 

204 

11  = 1 

242* 

K = 1 

243* 

IF (HO. LT. 01  PFTE=0. 

244« 

IFIMO.lt. 01  GO  TO  JOl 

24S>* 

ULE=ALOf,10(SUM(K  .MO.  I I .M.  IFLaGI  I 

246* 

THCOiiE:EXP(-16.fc*COH*UOI 

247* 

IRE=TRCONE«PCUnIC (ULE.NI 

24ft* 

0FD£=bETA(T0.V1-BETa(T211  1 .vmtre 

240* 

«R1 TE I6.60UIULE. TrE.OFHE 

2b0* 

OF1E:0. 

2‘jl* 

T2=T2(1) 

252* 

c 

computes  the  o<twi(  flux  for  the  layer  from  the 

253* 

c 

the 

reference  level 

254* 

IF IMO.NE .0)  GO  TO  JOS 

255* 

TRI I=TRE 

256* 

0B2=0B(T2. VI 

25/* 

CO  TO  302 

25ft* 

SOS 

00  300  HfI.MO 

259* 

I2=M«1 

2bU* 

K=K«4 

2bl  • 

T1=T2|KI 

262* 

IF Im.oT. 1 1 GO  TO  30G 

263* 

THIISTRE 

204* 

Oli.^FOnl  T2.VI 

2o5* 

30>> 

ULIFALOCIOISUMIK.MO. I2.M. IFLAGI  1 

r 

I 


49 


i 

i 


2bl>» 

ioO* 

eio* 

2Tl« 

27J. 

27m 

27!>. 

27.,. 


2 7'/» 
2b0« 
2B1» 
282* 
28J* 

2B‘»* 

2eb» 

2Bb« 

2b7. 

2A8* 

2e7« 

290* 

291* 

292* 

29J. 

29‘*« 

2Vb» 

29b. 

297. 

298. 

299. 
iuo* 
301* 

302. 

303. 
30-.. 
30b. 
30b. 
307. 
3U8. 
309. 
3l0» 
3ll. 
312* 
3l3. 
319. 
Sib. 
316* 
317. 
3ia. 

319. 

320. 

321. 

322. 


Ol:U0-U|K ) 

OBUDlUTl  .V) 

1»COf;i:'^*l’l-lb.6*CON»Un 
TH  IrTKCO'll  .PCUUIC  (ULl  .Nl 
DH  = .5«  IDri’.TPl  I •[•UMTRI  ) .(11-72) 

D7  lt=PF  lE.bFf 

.91  Te  (6.860  »K.  Tpl  .OFE.ULl  .U1  .0(11  »OFIE 

THI llTRI 

0B2IOB1 

300  T2:Tl 

302  IFLAGr3 

r COMPUTES  TMF  009IJ  FLU*  FOR  THE  La^ER  FROM  100  M ABOVE  THE  REFERENCE  LEVEi 
C TO  THE  REFEREbCE  LEVEL 
00  301  prl.S 
K:K.9 
Tl:T2(Kl 
«2rM.l 
U1FU0-U(KI 

307  IF (M2. LE. 51  GO  TO  309 
OLi:-3.8 
U1=0. 

GO  To  305 

309  UL1:alO610(SUM(k«mo.M>m2» IFLAG) ) 

305  0B1=DB(T1.V) 

TRC0m:E''F<*15*6*C0N.U1  ) 

THI:TrC0N1 .PCUblC(ULl.N) 

CFEf.5.(0B2»TRI l.nulATRI) .(71-72) 

OF It:OFIE*OFE 
ak2:k-4 

.HI T£ (6.660)k. Tri .OFE.ULl .Ul.OBl "OFlE 
IF  (K.lT.‘i69.0R.K.6T.>*77)G0  TO  309 
»R1Tei6.602)K.T2.CON.ThiI.OF1E.CFBE.0fe 
»RlTE(b,603)K,D01.0B2.0FlE 
«R1 TE (6.603)KK2. ri.ULl.TRI 
309  0B2rDBl 
TRI 1=TR1 

301  I2FI1 
0FTE:DFpE-0F1E 

c computes  the  up  flux  aT  the  reference  level  for  The  exact  case 

303  UFlErO. 

UFUEioETA(TO.V) 
tF(v,u,LT.O)  GO  to  903 

C COMPOTES  THE  UPFLUX  FOR  THE  LATER  FROM  THE  REFERENCE  LEVEL  TO  100  M bCLOW 
00  900  MFl.S 
K:K«9 
Tl=t2(X) 

U1=0(K)-U0 

IF (M.gT. 1 ) GO  To  909 

TRIln. 

T2sT2(I<-9) 

0H2=Dfl(T2.V) 

909  ULl:ALOr.lO(SUHU(K.MU.M.M.  (FLAG)  ) 

OH1FOH(T1>V) 

TRCOfilif  »P(-16.6*C0N.U1  I 
THUTHCPT'I  .PCHHlClULl  .N) 

UFEi  ,5*  (002.  TRI  1 .0.)l  .TRI  ) . ( TI-T2) 

UF  lEruFlf  .UFL 


I 

I 

i 


BEST  AVAIUSirCOPY 
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J?3. 

32b» 
3i?- 
32b* 
32'<* 
330  • 
331* 
33C* 
333* 
33‘«* 
33b* 
336* 
337* 
33a* 
339* 
3‘»u* 
3<4  1* 
3*42* 
3**3* 
3*4 ‘A* 
3'»‘j* 
34b* 
3*.7. 
34(1* 
349* 
3b0* 
3b  I* 
3b2* 
3b3* 
3b4* 
3bS* 
3b6* 
3b7* 
3b(>* 
359* 
3bO* 
3b  1* 
3o2* 
363* 
364* 
3bb* 
366* 
3b  7 • 
363* 
369* 
37u* 
371* 
372* 
373* 
374* 

37b* 

376* 

3/7* 

375* 

3/9* 


KK2:K-4 

44  1 Tf  (6. 660  IK  . If(I  .UFe>ULl  «0l  .001  •UFlE 
660  eOI.maT  I IX.  1 3.6eib.8) 

IF  (K  .L  7.4  7fl.0H.K  .r,  7 .40b(  GO  TO  407 

Mt  16.602  IK.  IH  I 1 .uEt  .gT  I E . I'ErE:  . 011  1 > u02 
»RI Tt (6.603IKK2. T 1 .ULl . TRI 
407  Tnll-THl 
002:01)1 

400  I2:ll 

If (MU.Eo.OI  GO  TO  403 

c computes  the  UPFLUX  for  the  LaTeP  FROM  100  M BELOW  THE  REFERENCE  LEVEL  To  TME 
C SURFACE 

IFLAG=4 

00  401  y:l.MU 
K:K*4 
TlrT2(K) 

406  ULl:ALOr,IO(SUMU(K.MU.M.y.lFLAGI  I 
Ul:U(K l-UO 
OBi:Oe(Tl .VI 

THCOl4l:EXP(-Ib.6*CON*Ul  I 
TRIaThCCNI *PCUOIC (ULl *NI 
4 0b  UFEf.5*  (0B2*rHn  *061*7017  * (T1-T2I 
UF  IE:UFIE*t'FE 

«RI  IE  (6.66n|K  . TRI  .UFE*ULl  .1)1  .OPl.UFIE 
TRI  I:TRI 
002:001 

401  T2:Tl 

403  UFrt:UFIE*HF0E 

402  *RITe16.601I 

601  FORMAT ( 1 X.6X. ’DfTC • . lOX . 'OcTE  • " lOX. 'OfTZ' . lOX, "UFTC* i 13X. "UFTE* . 

I lOX. .UFI2. I 

»R1 TE(6.60ni  0FTC.DFTE.0FT2.UFTc>UFTE.UFT2.0L0.V.ZRL.T0 
60J  FORMAT  ( 6 ( I X . £ 1 4 . 9 I . 4F 0 . 2 I 
GO  TO  502 
999  STOP 

FufJCTIOII  S>).((K2.MI  > 11 . J1 . IFLAC) 

SUM1:0. 

KS:k2 

IF ( IFLAG.E0.3I  GO  TO  6 
IF(Il.GT.Ml)  GO  To  10 
11  00  2 IS:  11. All 
00  3 JS:1.4 
KS:KS»l 

a2:10.**GL7(KS-14 

AIU0.*.GL7(X5) 

S:.5*(A2*All*  (UlKSI-U(KS-l) ) 

SUMli5UVl*5 
3 Sum: SUM I 
2 continue 

10  Ji:l 
6 004  IS:jl.S 
8 00  5 JS:1.4 
KS:KS*1 

A2:lO. ••GL7(KS-l I 

Alsi(l.**GL*(KSI 

S:.b* (A2*AI I • (U(KSI-U(K5-i 1 1 


BEST'AV/iLABLE'COPY 


jao« 

361* 

363* 
3a4« 
3ttb» 
30o» 
387* 
388  • 
369* 
390* 
39l« 
392* 
393» 
399  • 
39b» 
39b» 
397« 
396  • 
399* 
HOU* 
tOl* 
“Oi* 
403* 
404» 
405* 
40b« 
407* 
408* 
409* 
410» 

411. 

412. 

413. 

414. 
415* 
4|b. 
417. 
418* 
419. 
420* 
421* 

422. 

423. 
424* 
425. 
426* 

427. 

428. 

429. 

430. 

431. 
432» 
433. 
4 34. 
43b. 
43a. 


bOMlrSl'"!  «S 

IF(k2.LT.4<--,o.0H.K?.gT.477)  gO  TO  5 
»f(ITE  (6.bbniKb.Gl.2IKS-l  ) •GLZ(KGI  •42. 41.S>SU»'1 
650  F0H7'.AT(  IK.  I3.6(C15.8)  I 

5 SLI>.:51IF’1 
4 COflTIMUe 

9 RttORN 

FUNCTION  S""UtK3.M3. I 3. J3. IFLAG) 

SuM2rO. 

4U:K3 

IF  I IFL»r,.f  0.4)  GO  TO  4 
DO  2 I'J:l.J3 
00  3 jOrl .4 

Kgruu-l 

Ai:iO...GL7(KU»l) 

42=10. ••GLZIKUl 

Si.5.  (A2‘A)  I • (U(Ki|.l ) -UIKU)  ) 

SUM2:5Ur2.S 

lF(K3.LT.478.0R.K3.GT.4e5)  GO  TO  3 
»RIIE(6.650)KU,gL2(KU*1 ) . GL2 ( KU ) . 42 • A 1 . 5 . SUM2 
650  FOHWaTI IX. i3.6(e15.8I I 

4RI IE (6.655)U(Kg»l ) .glKU) 

655  FORI'^AK  lX,?E15.8l 

3 SUMU=SUk2 
IF(j3.E0.5I  SUM3=SURU 

2 COfiTIHUE 
return 

4 bUM2:SU'-3 

00  5 IU=1 . 13 

00  b JUz  1 . 4 
KUZKU-I 

41=10. ..CL'tKU*!) 

42=10. ..GLZlKUl 

S=.5. t A?.41 ) . IUIKU»1 )-UIKU) ) 

SUM2ISUV2.S 

6 SUKU=SU"2 
b CONTINUE 

RETURN 

FUNCTION  PCUBIC(X.,AR.NI 

1 = l 

IF(X0AR.LI.-3.7l  pCuOICrl. 

IF(*0AR.LT.-3.7I  return 
IF  (XRaR.GT . 1 .21  PcunlCrO. 

IFtXBAR.GT.1.21  return 
0X  = XBAR-XI 1 1 ) 

IF  lOXi  10. 30.20 

10  IF  (I.EO.ll  GO  To  30 
1 = 1-1 

0X:XMAR-X1 1 1 1 
IF  (0X1  10.30.30 

19  1:1.1 

Dx:0f)X 

20  IF  (Z.rO.NI  GO  To  30 
0DX:XI|AP-X’  I I . 1 I 

IF  lOOXI  3(1.19. 19 


^39* 
<«m  • 

(•44  • 
445* 
44b« 
44  7* 
Hh6* 
449* 
450* 
451* 
4b2* 
403* 
4544 
4bb* 
4b6* 
4S7* 
45fl* 
4i>9* 
4t>()» 
461  • 
4o2* 
4o3* 
4o4  • 
4o5* 
466* 
4674 
468* 
469* 
470* 
471* 
472* 
473* 
474* 
475* 
47to* 
4774 
4784 
4794 
4604 
481* 
402# 
4634 
4844 
46b4 
466* 

4674 

4664 

46‘>* 

490* 

49|4 

492* 

4934 


30  PBIC  =C(  I .1  )40X*(C<2.  I ) ♦0JC**C(3*P^0X*C(4i  I M ) 
PCUUIC=PBIC*.01 
RtTURN 

function  rmoczi 
IF  <Z-7?no.  .cr .0. ) GO  TO  I 

XX  = S4  22.Q9*  ( ( I ./2  73.  )-U  - /TDIZ)  ) ) ♦ALOG<fe»U  > 

E=te.XP(XX)  )MOOO. 

RhO=E/J  . 4614110. 4*7. I •T (Z» ) 

RETURN 

I XxX=bl47.50S4* ( ( I ./2  73. )-(  1./To(2»  M ♦ALOGCb.lU 
er(exp(xxxi  >41000, 

RmO=E/C  ,4G14U0.447.)4T(2)J 
RETURN 

Function  t<z> 

IF  (z-scio.  .r,E *0. 1 GO  to  i 
T=277. ♦(?. 75/500. )4Z 
RETURN 

I Tsaae. -.006542 

return 

Function  toiz) 

IF (Z-500..GE«0.)  GO  TO  X 
T0=274. ♦ (6./500. ) 42 

return 

1 IF  lZ-7200. .GE.O. ) 00  TO  2 
Tu=279.-<5fl.  /fa700)4«2-SOO.) 

RETURN 

2 T0=T(Z)-15. 

RETURN 

FufJCTION  P(ZI 

IF (Z-10769, ,gE.O. » GO  TO  I 

P=  C ( (Z^-44  3ne. > ♦! . ) 44 ( I ./. 190231 )4l0l3.25 

RETURN 

1 P=234 .52/EXP( (2-l0Tb9. >/638l.6> 
return 

FUNCTION  GL*TT#V*GL0) 

IF (V-1595. .GT.O. 1 GO  TO  I 
A:2. 4/(10. 445,) 

GO  TO  2 

1 A:l. 75/(10. 445, I 

2 6:2.22 

GL=GL0-A4( 1293, -TT)/TT) 4 C IV- 1595. )442. >4AL0GI 0(293. /TT) 
RE  ruRN 

Function  bftaitt*v) 

P0:3. 74) 1676/(10*4412. ) 

0:i .4389 

»ETAS(Pn*(V443,))/(EXP(04V/TT)-l.) 

R£  TURN 

Fuf.CT  IO»'  Ob(  TT.  V> 

Pit: 3. 741  1676/1  10.4412.) 

U:i.43R9 


BESI'AVAIWBIE'COPY 
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Oh:  (((Pf>*(TT*>2.))/tO*»3.))*(lO*V/Tr»**4 

49b» 

1 /TT  >-l 

<*9o» 

Ht  TURK 

497* 

49d» 

FUMCTIOH  07P(Tr .UZ.M.ZZpV.gLO.ZRL.CONI 

»99» 

500* 

IFIKZ.GT.O  ) GO  To  1 

bOl* 

JZ:ZZ-tO. 

iOi« 

IFtJZ.LT.O  1 GO  To  ? 

i03» 

IZ:ZZ-bon. 

bu4« 

IF  1 IZ.LT.  1 n.  AfiO.  IZ.GE.O)  GO  TO  2 

bi>3» 

IF  1 IZ.f.O.ANO.  IZ.OE.-IO)  GO  TO  1 

5Uo« 

221:22*10. 

b07. 

223:22-10. 

bus* 

TTl:T(221l 

bu9* 

TT2:TT 

bio* 

TT3:T(ZZ3I 

bn* 

IFIUZ.IE.O.I  OTRrO. 

bli* 

IFIUZ.LE.O.IGO  TO  3 

bli* 

UTI:gLITT1.V.gl0I *»LOGIO(UZ) 

bl4* 

UT3:gL<7T3.V.GLO) *ALOG10IUZ) 

bib* 

TCOM:EtP(-l6.6*COfl*UZ) 

blc* 

TRl:PCUnIC  tUT  1 .Ml  *10014 

bl7* 

TR3:PCUnIC(UT3.Ml •TCOM 

bis* 

OTR  :lTn3-TRll/(TT3-TTll 

bl9* 

GO  TO  3 

b2U* 

1 

TTUTT 

b2l* 

222:22-10. 

b22* 

TT2:T (2221 

b2J* 

IE(U2.LE.0. I OrR:0. 

b24* 

IFIUZ.LE.O.IGO  TO  3 

S2b* 

UTi:GLlTTl.V.OLOI*('.LOOin(UZI 

b2b* 

UI2:gL(TT2.U.GL0I *ALOGIO I U2» 

b27* 

TCOi::ExP(- 1ft.  6 *00*1*021 

b2s* 

TRlrPCUBlC lUT 1 .Ml •TCOM 

b29* 

TR2:PCUqIC(UT2.Ml *ICOM 

bJu* 

OTR  :(TR2-TRl I/(TT2-TT1 1 

bil* 

GO  TO  3 

b32* 

2 

TT3:TT 

bJJ* 

ZZ2:22*in. 

bJ4* 

TT2:T1227I 

bJb* 

IFlUZ.Lf.n.l  0TH:0. 

bJ6* 

IFIUZ.LE.O.IGO  TO  3 

bi7* 

UT2:GLITT2.v.gl0I *AL0Gin(u2) 

bie* 

UT3agL<TT3.V.GL0I*ALO&10(UZI 

bi9* 

TCOr.:EXP(-  Ift. ft *00*1*021 

b40* 

THZlPCOfilC  IUT2.M  1 • TCOM 

34  1* 

rR3:PcO«TC  (UT3>fi)  • TCOM 

b42* 

OTh  :(TP3-Tb2I/ITT3-TT2I 

b43* 

10  TO  3 

b44* 

3 

IF (ZZ.GT.bftO. .0R.Z2.lt .4ftO. IGO  TO  «0 

bHb* 

WRITE  Ib.bftOMT  1 . TI2.  TCOM.UTl  .(IT2.UT3.U2 

346* 

WRI  T£  (6.6G0I  tm  , TH2.TR3.  Til  .1  T2.  T r3.DTR 

34  7* 

660 

F0PFiaT(  IX.7(E  Ib.Rl  1 

b4S* 

4U 

HE  TuRM 

b49* 

ENO 

/( (EXf(0»v 
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.i'JunlobS6«Pl 

.5729b21-7*r,l 

. o'*  1645  7 3f.  *0 1 

• 7^J^J4o65^l9♦  n 1 
.750466535*01 
.a03565‘ir2*nl 
.659 170957*01 
.919731712*01 
.951902051  *r.l 
. 1062bb032*02 
. 1 1571*5726*02 
. 12b75U'0=!*0? 

•139250259*02 
. 153247737*02 
. lb875857S*ii2 
. 10571 5002*02 
.2o436952J*.',2 

.223rto3bP7*o? 

.2444i5u4o*c2 

•2b553372^♦^^ 

• 206*4lP,72*c2 
. JU7627C20*  ‘12 
.327770160*02 
.34j222353*r.2 
.3o334n32;*n2 
•3764101 14*02 
.390995u5:1*02 
.1*00603356*02 
.406591539*02 
.4u903o447*r2 

.406286624*02 

.432623000*0? 
.392421331*02 
.37649h,'>0*o? 
.354011655*02 
. J260606P6*02 
.291144607*02 
.252360659*02 
.211412759*02 
. 170980324 *„2 
. 133633972*0? 
,101475778*02 
. 754628990*01 
.556726.- 29*  01 

. ‘*05406160*01 

.2  »9'.  j9</05  *0  I 
.223o777‘-2*01 
. 15499991.5*  111 


LUi 


*^9  / 

r 


Interpolant  Coefficient  Input 


.49  10-13084  *01 
.51  734,1  756*0  1 
. 4 39536*;5  4 * 0 I 
.42451‘’371*oI 

.2o?4  1725'i*oi 

.325791  li^.O^Ol 
.23440833 1 *9 1 
.2’-i550  7'i  1*01 
.319274001 *ol 
.206332*458*01 
.333371003*01 
.472212577*01 
.475758100*01 
.h£477n558*ol 
.625l5«7‘  1 *01 

. 774560  >5  1 *0 1 

. 77t,H94ci3,,*ol 
.91 *444 i32*qi 

.9457o°01 9*0l 
.997720727*01 

. I CO336960  *02 
. 1043d79i7*02 

. 104 163003*02 

.100357991*02 

.623733756*01 

.0214.86294*01 
790316531 *01 
.717254060*01 
.54064--,  I 1 1 *01 
.420181298*01 
. I 78642^70*01 
.652503932*00 
.l39c?8J“4*ol 
.406731  1 .35*01 
.033433507*01 
.95952791,9*01 
. 122345697*02 
. 1626640.32*0? 
.186408618*02 
.201340742*0? 
.2o81382')6*o2 
• 196107.143*02 
. 177437301 *02 
. 144  l**i*o?8*l,2 
.1159.1.1100*02 
.8191  3*,iri  1*01 

.7.7,3  4141.7.1  *01 

.42715721 7*0) 
. 327  )5*24.,*01 
. nouo.liilio'1 


-.875228740*00 

. 259348 184*01 
.5007-1  947.4  *00 
.5403377,22*01 
-.21  1244717*01 
. 30'. 608378*01 
-. 72*477309-01 
-.275647670*01 
. 109808457*01 
-. 163459699*01 
-.4561 32936*01 
-.118! *6420*00 
-.4967047,32*01 
-. 129252884-01 

-.468049552*01 
-.o37c54509-01 
-.476479584*01 
-.8754731 19*U0 
-. 173387495*01 
-.2188  73048*01 
.460379240*00 
.22.4781001*00 
-.135802543*01 
.532817739*01 
.7495n5097-01 
.437230980*01 
.243549374*01 
.5»87o3042*D1 
.40I‘-44553*01 
.805149162*01 
,377079264*01 
.682645222*01 
.890350306*01 
.755682021*01 
. 105698840*02 
.8°54?5247*0l 
. 13272762-3*0? 
. 794461 765*01 
.494742393*01 
.226503472*01 
-.401 067  746*01 
-.622291535*01 
-.11 0675497  *0? 
-.9;8682520*01 
-.11 3551390*02 
-.519239879,01 
-.757522043*01 
-.  3.30672860*01 
. no?.*«4,'i5*0l 
.ocoooooon 


.993196997,02 
.980796  J9b*0? 
.9551.00000*02 
.980106002*02 
.9739069Ou,02 

.907296995.02 

.960000000*02 
.95?2c''003*02 
.9439o"003*02 
.930000000*02 
.925500002*02 
.915200997*02 
.904200001 *02 
.692X50000*02 
.873800001*02 
.864200001*02 
.843100004*02 
.830400000*02 
.810900002*02 
.769406998*02 
. 766096997*02 
.740606004*02 
.713006002*02 
.683299999*02 
.651496996*02 
.61 7800002*02 
.582300000*02 
.545200000*0? 
.50670000 1*02 
.467100000*02 
.426700001*02 
.585900002*0? 
.345000000*0? 
.304400001 *02 
.264600000*62 
.22h09699'7*02 
.189500000*02 
.155400000*02 
.124500000*02 
.973000002*01 
•741000003*01 
.550006000*01 
.3979<i6999*01 
.281000000*01 
.193006001*01 
.128000000*01 
.796Q4',o9  7«00 
.449900099*00 
. 18990699... 00 
.0000O9U00 


X. 

1 

-.36<5999«>99  + nl 
-.359999999*0l 
350000nr0^C  I 
-.340oonnoi*oi 
-.■i3nooonoi*oi 
-.31999«'^99*01 
“.30O999'»99  + yl 

-.30000nn00*0l 
-.290000npl*0l 
-.2sooonnoi ♦oi 

•.2699999y9^0l 

-.250999999  ♦•Ol 

-#250000000^0l 

-.2‘40000001*0l 

-.23GOOOrol^Ol 

-.2l®999999*nl 

-.2nO999999*0l 

-.2aoouoooo*’0i 

I900090ni^0l 

-.150000000*01 

-.i7nooonno*ci 

l5Q99Q«Q9*nl 

150000000*01 
140000001*01 
l3COOOroO*Ol 
-. 120000^00*01 
-.  10*5999099*01 
-.  lonooonoo^oi 
-.300999999*00 
-.799999007*00 

700000003*00 
-.630000001*00 
-.500000000*00 
-.399999'i99*00 
-.300000001 *00 
109990ouq400 
-.9999999M0-OI 
.000000000 
.990990054^-0  I 
, IQO999C99*.i0 
.300000001 *00 

.399999099*00 

.500000000*00 
.600000001 *00 
• 700000003*00 
.79O999n9/*00 

.090999099*00 
. 100000000*01 
.lO''999ooQ*Qi 

. I20000000*ul 


I 


Model  I Functions 


4bl* 

PUMCTICt;  T(2) 

1 

T:2'lfi.-.00<>*i*2 

HSh* 

4^  TUPf  1 

4t)6« 

FUflCTru^  TO(Z) 

*4i>7* 

1 

lF<2~7«‘'0n..rE.0.)  60  TO  2 

4b^* 

T6=2«3-  tS'i./bTOO. ) *2 

4b9* 

Ht  TUHn 

4bii« 

2 

4ul* 

KETUPM 

S • 

FV7I7.  T 1 1,*4  (itR(  1 1 r.^O.  ..'l  ,cn''l 

<<'J 

<<yt» 

If  (KZ.'iT.O  ) &;<  TO  1 

J7-22-10. 

49o» 

IF ( j2.lt. n ) GO  To  2 

9‘)7» 

221t,:£*U’. 

»9d* 

223:?2-’0. 

TT  ITTIZZl ) 

500* 

TIZrTF 

bul* 

TT3:T(7T3I 

bu«:» 

IFIUZ.ir.O.)  OTOzO. 

5J^» 

IFIUZ.LF.O.ICO  TO  3 

blj^  • 

UT  Iz.iL  IT''!  .V.GLO)  .ALOGIO  lUZ) 

5oS* 

UT3:c-i-<TT3.V.GLO)  .FLOGIOIUZ) 

bub» 

TCOI.rLXf'(-!6.6.CO.!.U2) 

b07- 

Tfil:PC'.'''ICHiri.i!)  .TCON 

50o» 

TR3FPCURtC  IUT3.III  .TCOM 

509* 

DTP  z(TR3-TRll/(IT3-TTl) 

blU» 

GO  To  3 

bll« 

1 

TTIZTT 

bl2* 

222=22-10. 

bl 

772=7 (222) 

Sli** 

IF(U2.LL.n.)  OTPrU. 

blS« 

IF (UZ.Lr .0. )G0  TO  3 

516* 

'JT  1 =gL  ( T T 1 , V . r,..0 1 . «L(TG  1 0 ( U2  1 

bl7. 

•jT2  = '.L<TTa.V.GLOI.ALOGlO(U2) 

5l.l« 

rc0li:FX'’(-|f..G.C0ri.iJ2) 

bl9* 

T;<i:pcL'°IC(UT1  .t.)  .TCO'I 

b20* 

7M2=PCUPrC('/T2<-|l.  I'COM 

b2l» 

3T,7  z(TP?-tpi)/(IT2-TTI) 

522* 

GO  TO  3 

S2J» 

2 

TT3=TI 

b2u> 

222=22.10. 

525* 

772=7(2221 

526* 

IF(u2.L€.n.l  UTRrO. 

527» 

IF(U2.LF.0.IGO  TO  3 

S25« 

ut2=gl(tt2.v.gloi  ..vLor.iocjz) 

b2'l» 

UT3=GL( TT3,v.GLC> .AlOGlOlUZI 

5J0* 

fC0/J=£»P(-16.6.C0iJ.i'2> 

bil« 

7K2rPClJP:C(ljT2.(il  .TCOM 

b32* 

Tr3=PCU0IC(UT3.;i)  .icon 

533. 

uTR  =ITR3-tr2)/|IT3-TT2) 

53'.. 

GO  TO  3 

935* 

3 

IF (22.GT.5G0. .On.ZZ.LT.u^O. IGO  70  40 

b36. 

WRI  7C  (6.660)77  1 . 7 T 2 . I COfI  > LiT  1 .l)T2  • UT  3 . U2 

537. 

2R7TC (6.660) 7hI . 7p2.7R3.T71. 772. TT3.0TR 

b33. 

660 

f 0»(.4aT(  i7,T(1  Ib.ftI  I 

b39. 

«0 

Rt  TURK 

5«u« 

zuc 

i 

1 
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